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Abstract 


The methanol oxidation reaction has been studied on polycrystalline platinum surfaces, cathodically treated platinum and platinum surfaces 
modified by a second (or a third) metal deposited spontaneously in 0.50M methanol +0.5M sulphuric acid using chronoamperometry. The 
electrochemical behaviour of methanol on the three types of electrodes shows that the initial oxidation rate is higher on the cathodically treated 
electrodes, whereas it depends on the nature of the second spontaneously deposited metal; silver shows inhibition and ruthenium and osmium 
exhibit a catalytic behaviour for the same 2096 amount of surface metal ad-atoms. 

Analysis of the methanol chronoamperometric data show that the steady-state current, recorded after 3000s, is larger for the cathodic treated 
surfaces than for bare platinum. Moreover, the initial slope of the current versus time plot is positive for both surfaces, that is, the competition between 
carbon monoxide formation and methanol first oxidation occurs. On the other hand, methanol oxidation on platinum modified by spontaneously 
deposited metals always exhibits negative signs in the initial current versus time slopes, independent of the nature of the metal ad-atom. 

A theoretical analysis of the kinetic laws for methanol electrocatalytic oxidation is presented. The fitting of the chronoamperometric data with the 
mathematical model, allows the discussion of methanol decomposition, carbon monoxide oxidation and direct methanol oxidation on the different 
surfaces. The time dependent expression of the surface coverage by methanol adsorbates is also derived and the characteristic values of the reaction 
are calculated for the different platinum electrodes. The potential dependences of parameters directly related to the electrochemical rate constants 


of the each step were also evaluated. 
© 2007 Elsevier B.V. All rights reserved. 
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1. Introduction 


Global warming due to greenhouse gases has been a world- 
wide problem since the late seventies. In order to decrease 
pollutant levels, advances in electric vehicles that produce 
negligible amounts of harmful emissions and show efficiency 
improvements are essential. Fuel cells are promising power gen- 
eration systems of high efficiency since they are able to convert 
the free energy change of a chemical reaction directly into elec- 
trical energy. Although an attractive strategy, it is obvious that 
improving the efficiency of current systems alone will not be 
sufficient to meet the Kyoto targets [1]. In the case of low tem- 
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perature fuel cells the selected fuel is hydrogen. However, the 
ability of hydrogen to compete economically as a transporta- 
tion fuel rests primarily on reduction of the fuel-cycle cost and 
storage methods [2]. 

One of the most popular liquid fuels is methanol, and there is 
great interest in its direct oxidation without a previous reforming 
step in the so-called direct methanol fuel cell (DMFC). One of 
the obstacles to overcome in making a DMFC practicalis the lack 
of efficient poison-tolerant anode catalysts [3]. In the oxidation 
of methanol on platinum, surface carbon monoxide is formed as 
a stable product, which blocks catalytic sites and inhibits fur- 
ther methanol oxidation. Therefore, a good anode catalyst for 
DMFC must transform the poisonous carbon monoxide into free 
carbon dioxide using oxygen atoms available in the interfacial 
water or in the water decomposition products such as adsorbed 
hydroxyl species. This leads to the release of surface sites for 
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a catalytic methanol oxidation. In this sense, platinum surfaces 
have been modified either with highly reductive potential pro- 
grams or with the deposition of ad-metals to change the carbon 
monoxide surface binding properties and lessen the effect of the 
carbon monoxide poisoning [4]. 

Although the decomposition and oxidation reactions of 
methanol on platinum have been studied extensively in the past, 
the exact mechanism that explains the observed experimental 
facts still remains unclear. It has been studied by Bagotzki and 
Vassiliev [5] on polycrystalline (pc) platinum at different elec- 
trode potentials, pH and bulk concentrations. They found that 
the decrease in the chronoamperometric transients were due to 
the increase in the surface coverage by adsorbed methanol. They 
also analysed the kinetic currents at different degrees of coverage 
by hydroxyl and carbonaceous species, making visible the sur- 
face recombination between the two adsorbates as rate-limiting 
step for potentials below 0.70 V versus RHE. Steady-state Tafel 
lines for methanol oxidation were used to analyse the change in 
the slopes considering different values of surface coverages by 
methanol. A similar analysis was performed by Indara et al. [6] 
for both methanol and hydroxyl adsorbates. 

Technological studies of methanol oxidation at different 
anodes in a DMFC shows that an increase in platinum load- 
ing to 5-10 mg cm7? is required to obtain the desired output 
current and potential in comparison with the 0.2-0.5 mg cm? 
on supported carbon catalysts [7]. The co-deposition of two 
or more active metals in a new binary alloy catalyst lowers 
metal loadings and then, the cost of the catalysts. According 
to previous papers [8-11] it seems that the most convenient 
electrocatalyst is the platinum/ruthenium alloy that can be pre- 
pared according to different methodologies; painting, spraying 
or ink printing, emulsions, electrodeposition, etc., over a car- 
bon support [12,13]. Stab and Urban [14] have patented an 
archive covering noble metal and binary alloy deposition meth- 
ods using sixteen metallic salts as anodes for a DMFC. It has 
been proposed that potential programs of different wave signals 
are the most adequate to obtain selected morphologies with fre- 
quencies down to 100 KHz. One simple method that does not 
require any sophisticated equipment is the spontaneous deposi- 
tion previously studied by us [15] and it is the one used in this 
paper. 

Considering the case of methanol oxidation mechanism on 
bimetallic surfaces, the “bifunctional mechanism" and the “lig- 
and effect theory" can explain the results. However, it has been 
demonstrated that the first mechanism contributes more than 
two-thirds of the overall effects found for methanol oxidation 
on platinum/ruthenium. Using the Polanyi relationship, it was 
possible to calculate that the process only changes in 40 mV due 
to the ligand effect of the ca. 200 mV shift in the onset oxida- 
tion potential towards the negative direction [6]. Gasteiger et 
al. [7] used different proportions of bulk platinum/ruthenium 
alloys in acid solutions to approach the chronoamperometric 
transients at low methanol concentrations. Other authors [12] use 
chronoamperometric plots for only qualitative purposes between 
the catalytic activity of bare platinum, platinum/ruthenium and 
platinum/tin. Herrero et al. [13] use the current at zero time 
as criteria in chronoamperometry to further obtain Tafel lines 


on Pt(1 1 1), Pt(1 1 0) and Pt(1 00) in sulphuric, perchloric and 
phosphoric media. They also analysed in UHV conditions with 
radiotracer methods the mechanism of methanol oxidation elec- 
trocatalysis with a kinetic law to model the chronoamperometric 
curves. The total current was considered as the addition of the 
double layer current contribution, continuous methanol oxi- 
dation component and seudocapacitive methanol adsorption 
current [12,13]. They solve the equations for different reaction 
orders for surface methanol species from one to four, founding 
the best fitting for the second order. 

On the other hand, a reaction mechanism has been proposed 
taking a single carbonaceous adsorbed species with further reac- 
tions including various types of hydroxyl adsorbates [16]. They 
assume that hydroxylates are fast and reversibly produced and 
then, their surface coverage can be almost constant at any reac- 
tion operation time. They fit their experimental curves with 
a proposed model finding good correlation. Ishikawa et al. 
[17] developed a theoretical study for methanol electrooxida- 
tion mechanism using the adsorption energy of the intermediate 
adsorbates by the bond energy of the metal atoms. They found 
that the first methanol dehydrogenation is favourable, but not 
the presence of water-discharged products, which are determin- 
ing in the oxidative desorption of carbon monoxide. They also 
found that the presence of adsorbed water on ruthenium atoms 
is more favourable than on platinum atoms. 

Based on the reaction scheme for methanol oxidation, it is 
clear that the ratio and distribution of platinum and ruthenium 
catalyst sites influences the methanol oxidation kinetics [18,19]. 
In spite of being the optimal surface distribution of both metals 
50:50 (as predicted by the Langmuir-Hinshelwood mechanism) 
it has been shown that catalysts consisting of homogenously dis- 
tributed surface sites of 70—30 atomic percent (at.%) platinum to 
ruthenium ratio display the best activity for methanol oxidation 
reaction [20]. 

The nature of the chemical state of the ruthenium ad-metal 
component has also been suggested to influence the methanol 
oxidation kinetics [21]; however, there is disagreement regard- 
ing the most active form of ruthenium. It has been reported to 
be significantly poorer (ca. 250 times lower) for methanol oxi- 
dation when catalysts consisted of mixtures of platinum and 
ruthenium-oxides [22]. Based on these results, it was recom- 
mended that platinum/ruthenium alloys should not be used for 
direct methanol fuel cells [23]. In subsequent studies [24,25], it 
was suggested that ruthenium in its metallic state is most active 
for methanol oxidation. They based their conclusions mainly on 
adsorbed carbon monoxide to carbon dioxide oxidation stud- 
ies for hydrogen-reduced catalysts. XPS was used to obtain 
the chemical state of ruthenium immediately before and after 
methanol oxidation [25]. It was shown that ruthenium-oxides are 
partially reduced to metallic ruthenium in the electrochemical 
environment [26] during the course of the reaction. 

However, it is important to notice that some earlier results 
on well-characterized platinum electrodes suggest that, under 
certain conditions, carbon dioxide production from methanol 
electrooxidation can occur by a non-carbon monoxide path- 
way i.e. the parallel pathway [27]. In a potentiostatic methanol 
oxidation it was also observed carbon dioxide production on 
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both Pt(1 0 0) and Pt(1 1 1) at potentials where adsorbed carbon 
monoxide is surface stable [28,29]. Direct evidence for carbon 
dioxide production at low potentials comes from in situ IR spec- 
troscopy measurements on both Pt(1 1 1) and Pt(1 00) electrodes 
[30]. The inferences from these studies were confined to poten- 
tials where adsorbed carbon monoxide did not oxidize to carbon 
dioxide. On the other hand, chronoamperometric measurements 
on Pt(1 11) at 0.6V versus RHE extended these results and 
showed that the parallel path was active even when adsorbed 
carbon monoxide oxidizes to carbon dioxide [31]. The parallel 
path accounted for more than 75% of the total carbon dioxide 
produced under steady state conditions [31]. In reality, methanol 
electrooxidation consists of several elementary reactions each 
of which may not only occur at different rates under unsteady 
state conditions, but also respond differently to residue accu- 
mulation. A model that successfully accounts for these transient 
effects must include further details of the surface reactions and 
not precisely using the parallel mechanism premise. 

In this paper, we develop a series kinetic model for methanol 
electrooxidation based on our current understanding. A the- 
oretical analysis based on the electrocatalytic kinetic law for 
methanol oxidation is presented. The fitting of the chronoam- 
perometric data is compared with the mathematical model and 
methanol decomposition and carbon monoxide oxidation are 
discussed on the basis of experimental and theoretical results. We 
analysed the model predictions with chronoamperometric data 
for acidic methanol oxidation on bare and cathodically treated 
platinum and on platinum modified by ruthenium, osmium and 
silver. 


2. Experimental 


Electrochemical experiments were performed at 20°C in a 
micro flux three-electrode cell using different pc platinum wires 
(0.5mm diameter, 99.999% purity from Goodfellow Co.) as 
working electrodes. The electrochemical set-up was completed 
using a large-area-platinum counter electrode and a reversible 
hydrogen reference electrode (RHE). The 0.5 M sulphuric acid 
supporting electrolyte was prepared from J.T. Baker reagent with 
MilliQ*-Millipore treated water (18.2 MQ cm of resistivity). All 
potentials in the text are referred to the RHE scale. 

Smooth platinum working electrodes were immersed first 
in hot concentrated sulphuric acid:nitric acid solution for half 
an hour, then repeatedly rinsed in Milli-Q” water and finally 
placed in the supporting electrolyte. Chemically cleaned plat- 
inum electrodes were electrochemical activated firstly by a 
20 min potential holding at 2.0 V in 0.5 M sulphuric acid solu- 
tion, secondly and after the solution was replaced by fresh 
supporting electrolyte, repetitive 2.5 V s^! scans were run for 
20 min and finally stabilized electrode by 15 min of potential 
cycling at 0.10 V s^! between 0.035 and 1.50 V. The real sur- 
face area of working electrodes was ca. 0.30 cm?, which were 
calculated from the hydrogen ad-atom voltammetric integration 
after double layer correction. 

The 0.5 mM ruthenium-containing solution was prepared 
from RuCl3-3H20 (Aldrich-Sigma), the 0.5 mM osmium- 
containing solution from OsCl5-xH5O (Aldrich-Sigma) and the 


0.1 mM silver-containing solution from silver sulphate (Merck) 
all in 0.2 M sulphuric acid solutions. They were used as a source 
of ruthenium, osmium and silver aqua-complex, prepared with- 
out local large salt concentrations to avoid precipitations and 
polymerisation. For spontaneous deposition, the immersion time 
(considered as the time in contact with the metal-containing 
solution) for ruthenium was 120 s and for osmium and silver was 
60 s under a nitrogen atmosphere. The calculation of the degree 
of surface coverage (Ome) by the foreign metal on pc platinum 
reaches 0.20—0.25 ML and was conducted by the integration of 
the voltammetric profile of the interface in the hydrogen adsorp- 
tion region after double layer correction according to Eq. (1). 
[9] 
ay, = Du En a) 
H 


Qt,and Qy being the hydrogen ad-atom charge density before 
and after metal deposition, respectively.Sometimes a solution 
composed by 0.5mM ruthenium-cation - 0.5 mM osmium- 
cation in 0.2M sulphuric acid was used to prepare the 
ternary platinum/ruthenium/osmium surfaces. The immersion 
time (faep) was 60 s. The stabilization of the spontaneous depo- 
sition on platinum was performed, after washing the electrode 
with supporting electrolyte, by not more than six potential cycles 
in the hydrogen ad-atom region at 0.10 V s^! until repetitive 
scans can be obtained. 

The cathodically treated platinum working electrodes were 
prepared by applying a potentiostatic signal (PS) in 0.5 M sul- 
phuric acid. Platinum was subjected to this PS at a preset 
potential, E, of —1.0 V during 5 min under a continuous bub- 
bling of the solution. The excess of molecular hydrogen at the 
electrode was removed by 15 min at 0.05 V under vigorous nitro- 
gen bubbling. In all cases the repetitive voltammetric response 
of the resulting surface was run from 0 to 1.55 V at 0.10 V s7! 
immediately afterwards to check the surface area of the working 
electrode. 

The 0.5 M methanol - 0.5 M sulphuric acid working elec- 
trolyte was prepared from pure methanol (Merck reagent) with 
MilliQ" -Millipore water. Methanol solution oxidation was stud- 
ied by chronoamperometric plots during 3000s in the working 
electrolyte with a sample frequency of 2 Hz, at different selected 
potentials on each platinum surface. 


3. Results 


3.1. Platinum modified surfaces by cathodic treatments: 
methanol bulk oxidation 


We prepared a freshly reduced platinum active surface by 
cathodic polarisation. Burke and Hurley [32] discussed the effect 
of platinum activation by thermal pre-treatment. He pointed out 
that, three active states of platinum are involved, which explained 
the three anodic peaks in the first positive anodic sweep in acid 
solution (0.75, 0.45 and 0.25 V). 

Fig. 1 shows the current versus potential response of a cathod- 
ically treated surface with PS of different E from —1.5 to 
0 V. The voltammetric response shows the following unusual 
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Fig. 1. Cyclic voltammetric profile of pc platinum (dotted line) run between 0.05 
and 1.55 V at 0.10 V s7! in oxygen free 0.5 M sulphuric acid solution at room 
temperature. The CT procedure is applied to the platinum electrode as explained 
in the text; E=0 V (continuous line), E- —0.5 V (dark grey line), £2 — 1.0 V 
(light grey line), E- —1.5 V (dashed line). 


features in the positive going potential sweep; a total deple- 
tion of the hydrogen ad-atom region and three anodic signals 
assigned to hydrogen and hydroxyl subsurface oxidation formed 
on more labile platinum active sites of lower coordination indices 
(located at ca. 0.3, 0.6, and 0.75 V) [33]. Moreover, at 1.2 V 
another anodic peak is observed associated with an oxygen- 
containing species of a less labile platinum site that does not 
allow the surface passivity. All these signals totally disappear in 
successive anodic potential scans, except for peaks at 0.75 and 
1.2 V that needs more than ten cycles to reach the usual repetitive 
platinum potentiodynamic profile. 

On the other hand, chronoamperometric plots at 0.70 V in 
0.50 M methanol + 0.5 M sulphuric acid on the pc and PS plat- 
inum is shown in Figs. 2 and 3, respectively. Methanol oxidation 
chronoamperometric decays were normalized with respect to 
the real surface area on PS treated platinum. The current densi- 
ties were larger than those obtained with pc platinum or metal 
modified platinum, denoting a strong catalytic behavior towards 
methanol oxidation as found before [54]. The open configura- 
tion of the resulting surface depicts the possibility of having 
new active sites for catalytic purposes [54]. This type of sur- 
face involves new active sites towards methanol oxidation not 
developed by metal modifiers. 


3.2. Platinum modified surfaces by a second metal: 
methanol bulk oxidation 


3.2.1. Platinum modified by inhibiting metals: silver 

The application of a potentiostatic or galvanostatic pulse 
is required in all electrodeposition methods; however, good 
results (especially at low coverages) have been obtained using 
the spontaneous deposition [34—38]. Large values of coverage 
can be reached by multiple spontaneous depositions as it has 
been demonstrated for ruthenium on Au(1 1 1) and for tin on 
pc platinum [39]. Besides ruthenium and osmium, there is lack 
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Fig. 2. Chronoamperometric plots for methanol oxidation on bare platinum 
recorded at different electrode potentials in oxygen free 0.50 M methanol + 0.5 M 
sulphuric acid; 0.40 V (continuous line), 0.65 V (light grey line), 0.70 V (grey 
line), 0.775 V (dotted line). 


of information about the physicochemical properties of sponta- 
neous deposition [40]. Silver deposition has been studied by a 
large number of techniques [41-44], being most of them on gold 
and platinum in sulphate or perchlorate electrolytes. Two differ- 
ent stages for silver deposition on platinum were observed; one 
at 1.1 V responding to silver—platinum alloy electrodissolution 
(overlapped with the oxygen-electroadsorption at free platinum 
sites) and the other at 0.65 V due to the silver oxidation (from 
the onset of the bulk deposition process) deposited on the former 
surface alloy [45]. 

The electrochemical contour of the silver spontaneous depo- 
sition can be observed in Fig. 4 from 0.05 to 1.50 V in oxygen 
free sulphuric acid solution. Spontaneous deposition produces 
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Fig. 3. Chronoamperometric plots for methanol oxidation on cathodically 
treated platinum (at E=—1V for 5 min in acid solution) recorded at differ- 
ent electrode potentials in oxygen free 0.50 M methanol + 0.5 M sulphuric acid; 
0.40 V (continuous line), 0.65 V (light grey line), 0.70 V (grey line), 0.775 V 
(dotted line). 
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Fig. 4. Cyclic voltammetric profiles of stabilized silver spontaneous deposition 
(dashed line) on pc platinum run from 0.05 to 1.50 V at 0.010 V s^! in oxy- 
gen free 0.5 M sulphuric acid solution at room temperature. Continuous line 
depicts the cyclic voltammetric profile of bare pc platinum. Silver spontaneous 
deposition is attained in 1 mM silver sulphate + 1 M sulphuric acid solution for 
tdep = 60s at room temperature. 


low surface coverages (Ag —0.20—0.30) comparing to those 
observed for the electrochemical deposition (not shown here). 
However, higher values can be obtained after successive spon- 
taneous depositions. On the other hand, the weakly adsorbed 
hydrogen peak is positively shifted from 0.12 to 0.02 V, whereas 
the "third anodic peak" (at 0.195 V) and the strong hydrogen 
adsorbate are similar populated by hydrogen ad-atoms. The 
main difference arising from higher coverage values of silver 
deposition is the appearance of a sharp peak beginning at 0.5 V. 


3.2.2. Platinum modified by catalytic metals: ruthenium 
and osmium 

Ruthenium and osmium depositions on platinum are of spe- 
cial interest to methanol oxidation in fuel cells. Among the 
different ways to deposit ruthenium or osmium on platinum 
[46—49], that of spontaneous deposition is attractive because of 
its simplicity and the fast surface concentration plateau reached 
after several seconds. They are very stable and normally change 
to stable hydroxides and oxides when the electrode potential is 
increased. The ex situ examination of the process by STM images 
on Pt(1 1 1) at the submonolayer levels shows the formation of 
islands [50] with maximum 0.20 ML coverage for ruthenium (at 
1205 of exposure) and 0.15 ML coverage for osmium (at 60 s of 
exposure). Depending on the osmium coverage the island diam- 
eter varies from 2 to 5 nm [51]. The island density increases with 
ruthenium coverage values, but in the case of osmium there is 
no optimum island size [51]. 

To check for the voltammetric features of foreign metals oxi- 
dation on platinum, voltammetric runs were performed with 
upper potentials larger than 0.80 V. Fig. 5 shows the cyclic 
voltammetric profile of ruthenium spontaneous deposition on pc 
platinum when the upper potential limit is extended up to 1.50 V. 
As found before, ruthenium species partially inhibits the hydro- 


L L L 1 fi 1 
0,0 0,2 0,4 0,6 0,8 1,0 1,2 1,4 
ENV vs. RHE 


Fig. 5. Cyclic voltammetric profile of spontaneous deposition of ruthenium 
(dashed line) performed from 0.05 to 1.50 V at 0.1 V s! (first and sixth cycles) 
compared to pc platinum (continuous line) in 0.5 M sulphuric acid. 


gen electrosorption reaction. Fig. 5 shows the slow increase of 
the oxidation current near 0.6 V probably due to the formation of 
RuO» and RuO; surface species detected by XPS measurements 
[52]. After extending the potential towards more positive val- 
ues, the partial inhibition of the oxygen electrosorption reaction 
is clear as a lower reduction charge density is observed in the 
reverse scan. Moreover, a sudden increase in the current previ- 
ous to the oxygen evolution reaction is also depicted in the same 
figure due to co-adsorbed chloride oxidation. 

Fig. 6 shows the adsorption of osmium-containing species on 
pc platinum when the potential limit is extended up to 1.5 V at 
0.10 V s-!. In the first positive going potential scan, the inhibi- 
tion of the oxygen-electrosorption is also evidenced, due to the 
presence of metallic osmium and osmium-oxygen containing 
species at a sub-monolayer level. According to results found by 
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Fig. 6. Cyclic voltammetric profile of spontaneous deposition of osmium 
(dashed line) performed from 0.05 to 1.50 V at 0.1 V s-! (two first cycles) 
compared to pc platinum (continuous line) in 0.5 M sulphuric acid. 
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XPS core-level spectra [53] and contrary to the case of ruthe- 
nium, metallic osmium can be found in the whole potential 
range, i.e., also up to 1.4 V. The nature of the oxidized species 
of osmium is OsO; in the entire potential region from 0.05 
to 1.50 V, but with only 25% of a metallic osmium monolayer 
[53]. Moreover, it is also seen in Fig. 6 that the oxidation peak 
of osmium species is only observed in the first positive going 
potential scan, probably due to the dissolution of Os?* ad-atoms 
at 1.0 V. During the second cycle the shoulder just prior to oxy- 
gen electrosorption reaction, located between 0.70 and 0.85, is 
better observed. 

The electrochemical response of ternary platinum/ruthenium/ 
osmium systems in the hydrogen ad-atom potential region was 
studied previous to methanol oxidation. The comparison of each 
voltammetric profiles run from 0.05 to 0.80 V at 0.10 V s! for 
the different surfaces shows that there are similarities between 
voltammetric profiles of surfaces covered by osmium containing 
species (not shown here). The broad pattern of the hydrogen 
potential region and the lack of well defined peaks for the three 
curves resemble that of a solely osmium deposition. It seems 
that osmium species governed the deposition process when it is 
deposited before or after ruthenium. The values of 0 resulting 
from the calculations of the voltammetric hydrogen ad-atom 
region were always ca. 0.25. 

We have performed the chronoamperometric curves in 
0.5M methanol+0.5M_ sulphuric acid solutions on bare 
platinum, platinum/silver, platinum/ruthenium/osmium, plat- 
inum/osmium and platinum/ruthenium surfaces (Fig. 7). The 
selected potential was 0.70 V, where oxygen-containing species 
(but no bulk platinum oxides) are present in all the studied 
surfaces. A background experiment was also performed in the 
absence of methanol (only with supporting electrolyte) to check 
forthe double layer-charging component at each potential. When 
subtracting this current contribution, we can calculate the early 
stages of the methanol oxidation rate at 0.70 V. 
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Fig. 7. Chronoamperometric plots for methanol oxidation on modified platinum 
(at E= —1 V for 5 min in acid solution) recorded at 0.70 V in oxygen free 0.50 M 
methanol 4 0.5 M sulphuric acid; bare platinum (black continuous line), plat- 
inum/silver (dotted line), platinum/ruthenium (light grey line), platinum/osmium 
(dashed line), platinum/ruthenium/osmium (dark grey line). 


It seems that the best catalytic effect can be attributed to bare 
platinum, but at very long times platinum/osmium surfaces seem 
to depict larger currents. The ratio between the slope of the cur- 
rent transient curve at zero time and the observed current at long 
times shows interesting results. The platinum/osmium surface 
exhibits the highest values of methanol oxidation rates, fol- 
lowed by platinum/ruthenium/osmium. The “synergetic effect” 
between ruthenium and osmium is evidenced in these experi- 
mental conditions. 

Theoretically, methanol dissociation on ruthenium (not in 
the case of tin or silver) is as favourable as on pure platinum 
as explained above. Water dissociation on ruthenium is sig- 
nificantly more favourable than on platinum. This means that 
we have the same probability for water decomposition on both 
surfaces but only platinum has the ability at room tempera- 
ture to adsorb methanol. The water dissociation to hydroxyl on 
ruthenium is slightly exothermic (—0.05 eV) in contrast to disso- 
ciation on pure platinum, which are quite endothermic (0.81 eV) 
[17]. Therefore, we expect that at the early stages of the process 
the bare platinum surface would exhibit a good performance 
upon methanol oxidation, but it has to decrease in comparison 
with the modified platinum surfaces; i.e. this is shown in Fig. 7 
where at times lower than 100 s, the best catalyst seems to be plat- 
inum, while at times ca. 800s the best performance is observed 
for platinum/osmium at 0.70 V as explained above. 

A totally different situation is observed for methanol oxida- 
tion on bare platinum surfaces at potentials larger than 0.7 V, 
where instead of having a current decreasing with time we 
Observe a current that increases with time (Fig. 2). This means 
that adsorbed carbon monoxide oxidation is faster than methanol 
decomposition reaction. The recorded current increases due to 
a slow increase in carbon monoxide coverage, which is rapidly 
oxidised. In this situation, the initial carbon monoxide oxida- 
tion reaction is fast but it slows down after some time as a result 
of an increased poisoning of the surface. Therefore, at poten- 
tials higher than 0.775 V, the incorporation of carbon monoxide 
oxidation reaction (not as poison but as a parallel reaction) is 
indebted, i.e. the potential range for our model is narrower due 
to this problem. Thus, an improved parallel mechanism together 
with an analytical solution of the differential equation is in 
advance. 


4. Methanol oxidation reaction mechanism 


In the platinum/aqueous acid solution interface and under 
0.4 V the hydrogen adsorption/desorption reversible reaction 
takes place: 


Ht +e7 © Hads (1) 


On the other hand, adsorption of water in the double layer 
with a further discharge from 0.7 V occurs in acid solution on a 
bare platinum surface: 


H20 + Pt <> Pt-H2Onas (II) 


Pt-H5O4q, <> Pt-OHads + Ht +e7 (IIL) 
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The interaction of these species with platinum is rather weak, 
except at potentials lower than 0.10 V where methanol adsorbed 
species are not able to displace hydrogen adsorbates and at 
potentials larger than 0.90 V methanol adsorption strongly com- 
petes with oxygen-containing species on platinum. 

The first interaction of methanol with platinum can be charac- 
terised by the subsequent losses of hydrogen from carbon atom, 
occupying various platinum adjacent sites. The number of these 
sites depends on the crystal orientation of the substrate, methanol 
concentration and pH. 


CH30H + 2Pt — Pt-CH5-OH,g; + Pt-Hads (IV) 
Pt-—CH2—OHads + 2Pt — Pt2z-CH—OHaas + Pt-Haas (V) 
Ptg=CH-OHags + 2Pt — Pt3=C—OHads + Pt-Haas (VD 


However, it has been shown that all of them are converted to 
Pt-CO (or Pt-COH) and some Pt-H depending on the involved 
potential range: 


Pt3=C—OHads > Pt=COaads + Pt—Hads (VII) 


The species Pt-COH and/or Pt-CO can react with surface 
hydroxyl species (III) when they are located at a nearby position. 
Finally, the products are oxidatively desorbed from the surface 
making platinum available for further acts of reaction. 


Pt3=C—OHads + Pt-OHads > Pt-—COads + H2Oads + 3Pt 


(VID 
Pt-COnds + Pt-OHaas > Pt-COOH,g, + Pt (IX) 
Pt-COOH,a, + Pt-OHags > CO2 -- H20 + 2Pt (X) 


A theoretical study has been performed using the exten- 
sive density-functional theory for the analysis of the adsorption 
behaviour of carbon monoxide on platinum and mixed Pt-M 
clusters (M 2 ruthenium, tin, germanium) [17]. This study exam- 
ined the effect of M alloying atoms in water dissociation, 
where the presence of M atoms decreases the PtM—CO bond 
strength. The activation energy for this process is smaller on 
ruthenium than on platinum. On the simulated cluster alloys 
of (Pt3)(RuzPts) water dissociation is also more endothermic 
and has slightly higher activation energy than on pure platinum 
[17]. On the other hand, methanol dissociation to COads occurs 
easily for pure platinum. They found that ruthenium in plat- 
inum/ruthenium alloys promotes methanol dissociation and the 
formation of hydroxyl adsorbates from adsorbed water. Their 
conclusion establishes that an alloy with a relatively high ruthe- 
nium/platinum ratio should be beneficial for methanol oxidation. 

Some earlier paper [55] assumed that methanol oxidation 
does not proceed through a series of consecutive reactions, 
but through a parallel path mechanism, forming a poisoning 
species through it. Thus, carbon dioxide may be formed by the 
direct oxidation of methanol or by the oxidation of adsorbed 
carbon monoxide at higher electrode potentials [13,28,56]. 
Using a chronoamperometric analysis of methanol oxidation 
on Pt(1 11), Pt(1 10), and Pt(100), Franaszczuk et al. [12] 
modelled the poisoning effect, with only the direct methanol 


oxidation and the formation of COads through a direct reac- 
tion pathway. The expression for the methanol-related oxidation 
current recorded in these transients’ contained four fitting 
parameters. However, it seems to be quite empirical because they 
assumed a rate constant of the site blocking process (decomposi- 
tion of methanol) with a single number of surface sites to convert 
methanol into carbon monoxide. 

Lu et al. [57] modified the former model to fit the experi- 
ments for methanol oxidation on pc platinum. They added an 
exponential current decay to the expression of the current tran- 
sients to account for the fast observed current decrease at times 
lower than 40 ms. This was ascribed to species formed during the 
initial oxidation, other than adsorbed carbon monoxide, report- 
ing a noticeable improvement in the fit. Vielstich and Xia [58] 
establish that the use of the initial current density in chronoam- 
perometric data as a fitting parameter is somewhat confusing. 
They also found that methanol does not oxidise completely to 
carbon dioxide because of carbon monoxide blockage, so a more 
complicated interpretation of the parameter has to be sketched. 

Stuve and co-workers [28] proposed the only model that 
incorporates carbon monoxide oxidation along with methanol 
direct oxidation and decomposition reactions. However, rather 
than using directly the chronoamperometric data, they chose 
to model the relationship between the charge under the tran- 
sients and the charge obtained after stripping the electrode from 
adsorbed species in a different cell. The results clearly showed 
that at short times (<30 ms) only partial oxidation products and 
virtually no carbon dioxide are formed. At times longer than 
5s, carbon dioxide production is greatly favoured over residue 
formation. They also concluded that for methanol oxidation on 
Pt(1 1 1) electrodes a serial reaction path, involving adsorbed 
intermediates, is inadequate to describe the observed rate of car- 
bon dioxide formation. Baltruschat et al. [59] proposed a parallel 
pathway mechanism in which carbon dioxide could be formed 
through the oxidation of adsorbed carbon monoxide or the oxi- 
dation of intermediate dissolvable species such as formic acid 
and formaldehyde. 


5. On the proposed model 


According to the above results by other authors, it seems that 
the parallel reaction pathway for methanol oxidation works for 
low plane and stepped platinum single crystals [60]. However, 
we propose for pc platinum a direct series oxidation mecha- 
nism setting the current produced in the direct path without 
making any assumptions regarding the nature of the reaction 
steps, the amount of electrons produced per methanol molecule, 
or the formed products. In other words, our model is sim- 
ple and without any improvements over previous models, but 
with a whole kinetic treatment for the overall series mech- 
anism. Nevertheless, many authors have shown [15,36] (for 
example) that the contribution of the direct pathway can be 
considerable on carbon monoxide-covered electrode surfaces 
and that the direct pathway has to be the main process of the 
mechanism. 

The following four-step reaction pathway in a series mech- 
anism is proposed neglecting the intermediates of methanol 


314 M. Ohanian, C.F. Zinola / Journal of Power Sources 168 (2007) 307—322 


dehydrogenation in the 0.4—0.8 V versus RHE potential ranges 
and considering the surface recombination of COads and OHads 
as the determining step. 


CH30H + Pr spp CH3_,OHags + xHt -xe^ XD 


Pt — CH3_,OHaas 


Foy zkia 2D L COpas + (4 — X)H* + (4 — xje^ (XI) 

Pt + Hosp; — Opa + H* + e7 (XIII) 
kox 

Pt — COggs + Pt — OHaas —> CO? + Ht +e + 2Pt (XIV) 


Kox,i and kreq,i being the electrochemical rate constants for each 
oxidation and reduction pathways, respectively. The rate deter- 
mining step, kox4, is of second order but since 0 is defined 
dimensionless, it is expressed in [s-1]. 

The following surface condition is applied: 


Oco + ĉo = 1 (2a) 
which is equivalent to: 


d6co _ oH 
t —— 8t 


In our model it is assumed that: 


(2b) 


(1) Thereis no distinction between each adsorbed carbonaceous 
species in the series mechanism. 

(2) All the oxygen-containing species are equivalent to 
hydroxyl adsorbates. 

(3) The reaction occurs under the Eqs. (XI)-(XIV) without any 
presence of platinum oxide precursors or hydrogen adsor- 
bates. 

(4) Bulk methanol and proton concentrations are considered as 
constants. 


5.1. The kinetic law of the series mechanism in methanol 
electrooxidation 


In the case of an electrocatalytic reaction; the current density 
(J) can be evaluated through the change in the surface coverage 
by carbon monoxide adsorbates (0) 


dé — dJ/Iu) Far dc 
d ^ d ^ NaAd d 


J = 0 (3a) 
c being the surface concentration in [mol cm ?], F the Faraday 
constant [C mol~!], o the maximum surface charge density in 
[C cm-?], I^ and I gat the surface excess at any time and at under 
saturation conditions in [molecules cm ?] with Isa =Nao/F, 
NA being the Avogadrois number. 

The current density for a Faradaic process is given by the 
electrochemical rate of reaction (v): 


JO = nFo(t) = nFK] [c^ (3b) 


n being the number of electrons, k — the potential dependent 
rate constant and p; the reaction order for bulk and surface 
concentrations of reactants. 

Considering Eqs. (3a) and (3b) and on the base of our model 
assumptions, we have the following expression of the kinetic 
law in the series mechanism: 


dé > j= 
aa; TEI — xkox,2 — Xkrea, (HT). — (4 — 3e o(HT) 7 
+ kred,3(H*)]0 + Fkox, 40^ — (4 — x) Fkox,ı(CH30H) 
+F Fkox,3 =0 (4) 
An interesting prediction that follows from Eq. (4) is that the 
shape of the current transient strongly depends on the relative 


numerical values of kox and keg. Therefore, by condensing the 
terms we will have: 


dé 2 
BER SR ak (5) 
with 


F 4— à 
= —[4- X)heao(H*) ^ — (4 — x)kox,2 + xkrea, (H^) 


— krea,3(H®)] (6a) 
L = Peas (6b) 
oO 
F 
M = [ko — (4 — x)kox,ı(CH30H)] (6c) 
lox 


According to (5) and (6) the dimensions of K, L and M are 
[cm? mol! s~!], whereas the dimensions of kox and kreq have 
to be [s^]. If the surface coverage were defined as an abso- 
lute coverage (T in [mol cm" ?]) the dimensions of k; would be 
dependent of the reaction order. 

Expression (5) is the quadratic Bernoulli equation, known 
as Ricatti equation, of non-linear differential equations with 
potential dependent coefficients [61]. Since the quadratic inte- 
gration of this equation is not analytically simple, the search of 
the solution cannot be reduced to a finite number of successive 
integrations. In Appendix A we have a brief summary of what 
we can do to solve the Ricatti's expression with the substitution 
Theorem. 

We can change to a new variable in Eq. (5) of the form: 

/ 
pep eA (7) 
ALT sat 

Then, Ricatti’s expression is reduced to a second order linear 

equation, which after reordering will be: 


Pat LM 
h” + Kh + = psn (8) 
NA 


The coefficients of the equation contain the potential depen- 
dent rate constants of the chronoamperometric experiments. 
Then, we can work with the constant coefficient characteristic 
equation. 

Tsai LM 


2 
Kr + —— =0 9 
r° + Kr + Na (9) 
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which has the following solution: 


K KM TLM do) 
T|,2 = zu 
1,2 2 2 NA 


The general solution of the differential equation will be: 

h = Cı exp(rit) + C2 exp(rot) (11) 
Taking the first derivative 

h/ = Cyr, exp(rit) + Cor? exp(rot) (12) 
and substituting it by Eq. (7), we obtain the value of 0 


|. WNa _ Na(Ciri exp(rit) + Cor; exp(rot)) 
ALT sat Isa L(C1 exp(rit) + C2 exp(r2t)) 


(13) 


and sinceL = FXox,4/o the current density expression will be: 


Wao! Nac? CiC»(ri — r3)! exp((r1 + r2)n) 
— 0 e 
J dt TyatFKox,4 (C1 exp(rit) + C2 exp(r2t))* 
24 exp((ri + r2)t) ; (14a) 
(C; exp(rit) + C2 exp(rat)) 


oC; C»(ri — r2}? 
Kox,4 


with Az 


since Isat = Nao/F (14b) 


This is the expression of the chronoamperometric current 
based on the series oxidation reaction model. Table | shows the 


C103 


r103 


0,40 0,45 0,50 0,56 0,60 0,65 0,70 0,75 
E/V vs. RHE 


Table 1 

Parameters for Eqs. (14a) and (14b) for chronoamperometric data at 0.70 V 
at different platinum and metal modified surfaces; ruthenium, osmium, ruthe- 
nium:osmium and silver in 0.50 M methanol + 0.5 M sulphuric acid 


A(Acm?) mnm m CA C2 
Pt 3.016 x 1076 — 0.06578 0.0527 0.12197  —0.0299 
Pt PS 5.606 x 1077 — 0.06977 0.05657 0.03254 — —0.00816 
Pt Ru 2.346 x 1076 — 0.06668 0.0348 0.13743  —0.04842 
Pt Os 1.533 x 1076 — 0.06922 0.0618 0.16037 —0.04111 
Pt Ru Os 2.085 x 1076 — 0.06958 0.0262 0.13758 —0.04723 
Pt Ag 2.958 x 107° 0.11816 0.0132 0.21208 —0.02542 


values of each parameter for the different platinum modified 
electrodes. 

The values of the fitting parameters at different electrode 
potentials were calculated for pc and PS platinum and plotted 
in Figs. 8 and 9. The values of the pre-exponential factor, A, 
includes the exchange current density of the reaction on both 
pc and PS platinum. It increases with the electrode potential, 
showing that the reaction requires the presence of hydroxyl- 
species on the surface to be activated from 0.6 V. On the other 
hand, the values of Cı changes with the potential in an oppo- 
site way as those of C2 as a consequence of the resulting sign. 
Cı is positive and has the lowest value at the largest studied 
electrode potentials, whereas C» exhibits its largest values (less 
negative) in the same region. This means that C, is involved in 
the methanol adsorption process, which is large in the hydrogen 
adsorption region (at potential higher than 0.10 V versus RHE). 


85r A 

80} (b) ^ 

75} a 

70} i 
65} 
60} EL. 
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Fig. 8. (a-d) Chronoamperometric fitting parameters: C2, Cj, r2 and rı for methanol oxidation on cathodically treated (open triangles) and bare (full squares) 
platinum. The cathodically treated platinum was prepared as explained in the text. The chronoamperometric curves were recorded at different potentials in oxygen 


free 0.50 M methanol + 0.5 M sulphuric acid. 
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Fig. 9. Chronoamperometric fitting parameter A for methanol oxidation on 
cathodically treated (open triangles) and bare (full squares) platinum. The 
cathodically treated platinum was prepared as explained in the text. 


C2 is concerned on the methanol oxidation process because of 
its potential dependence. Similar tendencies are observed for 
the PS platinum surfaces but with much larger values for C; and 
lower values for C5. 

However, one of the most important features in the cur- 
rent density expression is the shape of the current decay. This 
response can be envisaged through the values of rı and r». 
In the case of r2 the values softly decreases until 0.65 V, 
showing that this parameter is well associated to methanol 
oxidation. The values on PS platinum are larger until this 
potential value, but later are practically the same to those of 
pc platinum. On the other hand, the situation is not so clear 
for rı since the increase with potential seems to be saturated 
up to 0.60 V, coincidently with the minimum adsorption of 
methanol (competition with oxygen from the water discharge). 
The values are much higher for the PS surfaces at all poten- 
tials, showing that the current exponential decay is faster on pc 
platinum. 

A comparison between the fitting results in the case of 
methanol oxidation shows good correlation at times lower than 
200s, however at larger times the problem of overestimation 
appeared (not shown). One of the main reasons probably arises 
from the lack of correlation between our model and the carbon 
monoxide path in parallel. 

The complete expression of the current can be analysed 
through its slope at zero time, the current value at the asymptot- 
ically infinite time and the expression of the time dependence of 
the surface coverage by methanol adsorbates. 


5.1.1. The slope of the current density at zero time for 
methanol electrooxidation 

The slope can be evaluated from the current versus time first 
derivate: 


dj(t) E oCiCa(ri — r2 \ (ri  r2(C1exp(rit) + Ca exp(rat)) exp(Gri + r2)t) — 2exp((ri + r2)t) 


By taking the value at zero time: 


d(D0| | ([eCiComi = r2 
dt |,<0 kox,a(C1 + C2)’ 
2(Ciri + Cor2) 
x |(r1- r2) (16) 
| uM ee) 


The sign of the slope seems to be always negative in the 
base of this model, since the results of the fitting from Table 1 
show that (r1 * r2) &: 10(Cir; + C2r2) except for the only two 
exceptions of pc platinum with 5.5 1075 o/kox,4 and PS platinum 
with 6.9 1076 o/kox,4. These exceptions are more notorious at 
higher potentials as it can be seen in Fig. 3. 

According to the chronoamperometric data by Housmans 
and Koper [60] the transients recorded on Pt(11 1) showed 
an initial current increase, which may be explained by the 
carbon monoxide oxidation being faster than the methanol 
decomposition. This low decomposition rate is probably the 
result of a sufficiently low defect density on the metal and 
the low methanol concentration used in their experiments. 
They fit the chronoamperometric data with a mathematical 
model, which includes methanol decomposition, carbon monox- 
ide oxidation and direct methanol oxidation. The model indeed 
predicts that when the carbon monoxide oxidation rate is 
faster than the decomposition rate, a rising current transient 
can be expected, as we see for pure pc platinum. In this 
case it seems to be better to mimic methanol oxidation reac- 
tion with a parallel carbon monoxide and carbon dioxide 
mechanism. 

The analysis of the chronoamperometric parameters of 
Table 1 for metal-modified platinum surfaces shows that the most 
dramatic initial current versus time slope occurs on the ternary 
surface compound of ruthenium/osmium (—1.3 x 10-3 Olkoy 4). 
On the other hand, the lowest initial current versus time 
slope is observed for the platinum/silver surface compound 
with —7.4 x 1074 Olkox4. Besides, on platinum/ruthenium 
and platinum/osmium the slopes are —8.1 x 107^ o/kox,4 and 
—24 x 107^ Gkox 4, respectively. 

The analysis of the slope of the current versus time at 
t=0 involves the electrochemical rate constant of the deter- 
mining step (kox4) or said in other words the co-adsorbate 
catalytic oxidation, i.e. the oxidation of COads with OHads. 
Thus, depending on the value of kox, each slope can change 
its order of magnitude. In the case of having a large L, the 
values of the slopes will decrease. This means that the cur- 
rent observed for platinum denotes a low residential time for 
adsorbed carbon monoxide. On the other hand, a totally different 
situation arises for small L values; thus, the surfaces will denote 
a large residential time for COads, i.e. a slow oxidation of the 
adsorbate. 


dt Kox,4 


(C1 exp(rit) + C2 exp(rat))? 


(15) 
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5.1.2. The current density value at the asymptotically 
infinite time for methanol oxidation 

It is really important to evaluate the time at which the cur- 
rent transient reach the steady state, since it is a characteristic 
of the electrocatalyst. This is achieved considering in Eq. (15) 
the following reorder of terms by multiplying and dividing by 
exp(—(ri * r2)t): 


dj) _ (mem = rni + ro) expr + 2j 


dt kox,4(C1 exp(rit) + C2 exp(rat))? 
x (C; exp(rit) + Co exp(rot) — 2) (17) 
There is no change in 
dj(r) 
18 
di (18a) 


for platinum/silver is reached at the lowest t; values, show- 
ing again that the surface exhibits inhibitory characteristics 
(tss ca. 19 cte). 

(2) The other condition (18b) r; = r? implies 


K KM Peat LM 
nN =F] > 
2 2 NA 


K K\? DhaLM 
= + (22) 
2 2 NA 
This situation only accounts when: 
K\* DOLM 
>] = 23 
( 2 ) Na e 
that is, 


4—x 2 
NA F(xkrea, (H*) + (4 — Xx)kgeao(H*) — krea,3(H*) — (4 — x)Kox,2) 
= —AF? l'sa(x(CH3 OH)kox 4kox,1 — Kox,4kox,3) Q4) 


when either 


C; exp(rit) + C2 exp(rot) - 2 = 0 (18b) 


and/or r, =r 


(1) In the first case the comparison between both exponential 
terms can be reduced since they are both of the same shape 
but with a different expansion with time. Thus, they can 
be simplified to a single exponential contribution. However, 
looking at Table 1 we can say that rı = r2 +a, C1 =—(C2 +d), 
a and b being experimental parameters dependent on the 
type of platinum surface. Thus, Eq. (18a) renders for the 
steady state time, tss: 


C, exp(rifss) + C2 exp(ratss) — 2 
= —(C2 + b)exp((ro + a)tss) + C2 exp(ratss) - 2 = 0 


(19) 
or 


[C2 + C1 exp((r1 — r2)tss)] exp(r2tss) = 2 (20) 


The situation implies large times so we can say 
C2 « (C5 + b)exp(([ri — r2)tss), therefore: 


In(2 
orc Eis oe QD 
ri 


The relaxation time (at which the steady state is reached) 
indicates how fast the coverage of the intermediate will 
relax back to an equilibrium value after perturbation. Table 2 
shows the value proportional to tss obtained from Eq. (21). 
It is possible to conclude that for PS treated platinum the 
tss value is slightly larger than for pc platinum and more 
properly is the largest of all, i.e. ca. 59 + cte. On the other 
hand, the value of tss for platinum/ruthenium surfaces is 
near that of pc platinum, whereas it decreases for osmium 
and ruthenium/osmium surfaces. The steady state current 


However, this is a complex combination of electrochemi- 
cal rate constants. It would be better to analyse the difference 
between each rı and r2. On one hand, we can obtain the values 
of K, L and M from the difference between rı and r2: 


r-nV (KM Toa LM eS 
2 Et. NA 
On the other hand, from the addition of rı and r2, we have 
the value of K (Table 2) 


K = —(n 4 r2) (26) 


Considering that Eq. (26) renders negative values of K and 
since K is given by (6a) it is likely that the contribution of kox,2 
and krea,3 has to be larger than kyea,1 and krea,2. By substitution 
of K into (25) we can obtain the product LM. The values of L are 
given above from the slope of the current transient at zero time. 


ri- rV a qe +r) ^  DaLM (27) 
2 2 Na 
Then, by simplifying the expression we can obtain 
DhaLM (=r +r)? n-nWV Q8) 
Na 2 5 
or 


DaLM _ (rir) (29) 
Na 2 


The term r} + y calculated from Table 1 is always posi- 
tive and considering Eqs. (29) and (6) the product LM will be 
negative: 


cLM = —F^[(4 — x)(CH3OH)kox,4kox,1 — Kox,4kox,3] — (30) 


Since LM is negative, Eq. (30) indicates that ko 5 is smaller 
than (CH30OH)kox,1. However, we can separately obtain L from 
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Table 2 


Parameters K, L, M, kox 4, tss» Jt=0 and dj/dt, =o from chronoamperometric data at 0.70 V at different platinum and metal modified surfaces in 0.50 M methanol + 0.5 M 


sulphuric acid using data of Table 1 


jr=0 (A cm?) dj/dt, =o K LM LIF kox4 (s) | Mo (Amol!) tss (S) — cte 
(A cm ? s7!) (cm? mol-! s7!) (cmt mol? s72) (c7! cm? s7!) 
Pt 5.90 x 1074 7.68 x 10-6 —0.11848 —0.00155 0.2068 5 x 107? —0.00749 42.52 
Pt PS 1.00 x 10? 2.08 x 1075 —0.12634 —0.00167 0.0824 2 x 1075 —0.00749 59.05 
Pt Ru 6.62 x 1074 —1.97 x 1075 —0.10148 —0.00324 2.888 7 x 1074 —0.00112 40.16 
Pt Os 2.23 x 1074 —1.35 x 1076 —0.13102 —9.722 x 1074 0.563 1x 1074 —0.00173 36.46 
Pt Ru Os 5.27 x 1074 —2.26 x 1075 —0.09578 —0.00415 5.861 1x 1073 —7.0891 x 1074 38.47 
Pt Ag 2.24 x 1076 —1.13 x 1072 —0.13136 —0.01379 0.02 5x 10-9 —0.68938 18.99 


the initial current versus time slope (Eq. (16)), 
FCC — r2 
(Ci + C2) (dj(0/dr)io 
2(Ciri + Cor2) 


x 171 +r) (31a) 
Ov €i CD 
or in a simpler way from Eq. (14b) 
FCC — r2 
iz fel ari = r2) (31b) 


A 


Summing up, it is possible to obtain the values of L, M and 
K, which are depicted in Table 2, calculated from Eqs. (26), (29) 
and (31b). The most interesting parameter is L since it simply 
involves the value of kox,4. The largest value is observed for 
platinum/ruthenium/osmium surface (L=5.861F), so since L = 
Fkox 4/0. then, kox4 75.8610. Besides o 22.4 x 1074 C cm~? 
[60], thus, Ko 4 2 0.0014 mol! cm? s71. 

Table 2 also shows the values of ko, 4 for rest of the surfaces. 
At this point, it is important to remark that the values found for 
this oxidation rate constant by other authors [60] lie within the 
range of 10? to107! s~!, considering it as a first-order reaction 
constant. This situation can be also envisaged from Fig. 7 (dark 
grey lines) for times larger than 900 s, where the chronoampero- 
metric plot exhibits larger currents on this surface. Surprisingly, 
the value of L for PS platinum is almost three-fold lower than pc 
platinum. The situation is quite similar for platinum/ruthenium. 
This means that this type of surfaces are able to react better 
at starting times but later gets more poisoned than the original 
pc platinum. The situation for platinum/silver is the worst as 
expected. 

The values of K for the six types of platinum surfaces are 
nearly the same (the lowest for platinum/silver). The parameter 
contains the first back reductive methanol adsorbate, the back 
formation of COads subtracted the back formation of OHads and 
the decomposition ofthe first partially dehydrogenated methanol 
conversion to COads. The sign on all surfaces show that the two 
latter’s are the largest on all the studied platinum surfaces. This 
means that both processes are not very distinct on each modified 
platinum surfaces and also not rate determining. 

On the other hand, the value of M involves the forma- 
tion of OHads species minus the first oxidative dissociation 
of methanol. In all surfaces the sign is negative showing that 
the formation of OHads is not very fast. This is an impor- 
tant comparison since the problem of OHads formation is later 


determining the oxidation of COads. Thus, it is likely that the 
electrochemical rate for the first oxidation stage would be slower 
than expected, however it is not rate determining as it was sup- 
posed in our model. Some interesting features have to be pointed 
out, Mo = —0.0075 for both pc and PS platinum. This means that 
both electrochemical reactions: kox,; and kox,3 does not change 
with the cathodisation procedure, i.e. the number of active sites 
for both methanol and OH species are the same. However, in 
the case of platinum/silver Mo = —0.689 and is the largest of all. 
This means that the formation of OHads and CH3_,OHads are 
not the fastest of all, but the methanol oxidation reaction tends to 
be limited there, or in its COads stable species. Between ruthe- 
nium and osmium, it seems that the formation of OHads is slower 
on platinum/osmium and maybe also the formation of the first 
adsorbed species (CH3_,OHads). 

Another interesting correlation is the potential dependence 
of K, L/F and Mo calculated from current versus time fitting 
parameters: C5, C, r2 and rı as above. The results for methanol 
oxidation on pc and PS platinum surfaces are qualitatively sim- 
ilar, so we are going to explain the results only for methanol 
oxidation on pc platinum. The values of K increase exponen- 
tially with the electrode potential (not shown) as expected for 
the decomposition ofthe first partially dehydrogenated methanol 
conversion to COads. Though, from 0.65 V it sharply falls down 
to the initial values denoting the influence of the OHads back 
formation. The same situation is observed for L/F, since after its 
maximum at 0.55 V it rapidly decreases attaining the saturation 
at L=0.08F from 0.7 V. This means that the electrochemical 
rate constant of the determining step reaches its optimum value 
at 0.55 V. On the other hand, Mo values are nearly constant 
(between —4 and —5 x 107?) for potentials between 0.4 and 
0.65 V, however from 0.7 V it sharply decreases to ca. —0.022. 
At these higher potentials Mo is still negative indicating that 
the formation of OHads is not very fast comparing to methanol 
dissociation. 


5.1.3. The time dependence surface coverage by methanol 
adsorbates for methanol oxidation (integration of the 
chronoamperometric curve) 

On the other hand, by integrating the time dependent cur- 
rent density we can obtain the surface coverage by methanol 
adsorbates. Thus, from Eq. (14a): 


dð ^ Nae  CiCo(ri — r2) exp(ri + rdt) 
dt Tat FRox,4 (Cy exp(rit) + C2 exp(r2t))” 


(32) 
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[ ap = Cim = Y ri exp (ri +12) 1) p 
0; Kox,4 o (Cj exp(rit) + Co exp(rot)) 
(3 


or we can factorise the exponential terms to reduce their expo- 
nential order: 


[ qo CO = 0X [ 
Oi 0 


kox,4 


exp(rıt) 
(C2 + C1 exp(ri — r2) 


By defining 
Ciri — 12)" 

gan 21 — r2) (35) 

Kox,4 

we can obtain 
t 
t 
6 = 6; — % n a (36) 
o (C2 + Cı exp(ri — r2)t) 


0j being the initial surface coverage. Then, after multiplying and 
dividing by exp(ri — r2)t and exp— (rı +72)t 
t 
= t 
o (Ci + Coexp(ri — r2)t) 


We can use the substitution Theorem for integration as in 
Appendix B to obtain the time dependence of 0 


Ci(r1 — r2) | 1 1 


(37) 


0 — 6; 


(Cj; +C2) (Cı + Coexp(ri — r2)t) 
(38) 


kox,4 


The surface coverage at the beginning of the current transient 
will be as expected: 


lim |,.,00 = 6; (39) 
and the value of 0 at steady state (f — 00, Oss) is: 


Ci(ri — r2) 

lim |0 = Oss = 6j kon a(Ci + C2) (40) 
We can compare this magnitude defined as A0 = 0, — 0;—0 

in Table 3 using data of Table 1 taken for bulk methanol oxi- 

dation at 0.70 V on the different catalysts. Since the change in 

0 strongly depends on koX 4 we have to take into account that 

their values from Table 2 to calculate A0. The largest of A0 is 


Table 3 

Values of A0ko, 4 (with A0 = 6-9 — Oss), AO and Iss (at 1800 s) from parameters 
of Tables 1 and 2 for chronoamperometric data at 0.70 V at platinum and metal 
modified surfaces; ruthenium, osmium, ruthenium:osmium and silver in 0.50 M 
methanol + 0.5 M sulphuric acid 


AO0 kox4 (S71) A0 Js, (A cm?) 
Pt 0.0173 1.58 x 1078 3.1 x 1071! 
Pt PS 0.0176 3.91 x 1078 7.0 x 10-H 
Pt Ru 0.0492 1.68 x 1076 3.0 x 102! 
Pt Os 0.0099 3.82 x 10? 8.1 x 107? 
Pt Ru Os 0.0061 6.50 x 10-5 2.7 x 10727 
Pt Ag 0.1200 1.28 x 10? 1.3 x 10756 


observed for platinum/silver surfaces, meaning that Oss is very 
different to 0;..9. This does not mean imply large kinetics but 
a strong poisoning effect, demonstrated by the value of the sta- 
tionary current at 1800s, Zs. However, the problem is that 0; 
for this surface can almost be neglected and the change in the 
surface coverage covers insignificant values. On the other hand, 
in the case of PS platinum the value A0 is more than three-times 
larger than the reaction on pc platinum and also the value of Iss. 
This means that the consumption of adsorbed species during the 
process on these surfaces is very important, however the value 
of kox 4 is lower (also three times lower) than on pc platinum, 
so a compensation effect is expected. Summing up, the product 
AOkox 4 on both platinum surfaces is nearly the same and also 
the accumulation of adsorbed species on platinum/ruthenium 
and platinum/osmium, as it has been expected from the simi- 
lar coverages obtained by each metal ad-atoms. To make further 
conclusions on the comparative value of 0 a more advance model 
with a distinction between the initial methanol adsorbate and that 
of poisoning species is required. 


6. Conclusions 


(1) Chronoamperometric data were compared to envisage the 
catalytic effects of the different platinum surfaces towards 
methanol oxidation. The largest initial current density was 
observed for the cathodically treated surface; however, 
the asymptotic current showed a better response for plat- 
inum/osmium and ternary platinum/ruthenium/osmium. 

(2) An initial current increase at 0.70 V was observed in the 
chronoamperometry for pc and cathodically treated plat- 
inum, because of the competitive initial carbon monoxide 
oxidation and methanol decomposition. However, plat- 
inum surfaces modified by a second metal spontaneously 
deposited always showed current transients with negative 
slopes, independent of the nature of the second metal. This 
suggests a possibility for minimising carbon monoxide for- 
mation. 

(3) A theoretical analysis of the kinetic law for the proposed 
series mechanism of methanol electrocatalytic oxidation is 
presented here with an analytical solution of the differen- 
tial equation solved using Ricatti's equation. The solution 
obeys a linear combination of exponentials (one of them 
related to carbon monoxide formation and the other to its 
oxidation). The fitting of the chronoamperometric data with 
the mathematical model, leads to an analysis of methanol 
decomposition and carbon monoxide oxidation. 

(4) The time dependent expression of the surface coverage by 
methanol adsorbates is derived from the chronoampero- 
metric data. The asymptotic value of the surface coverage, 
A0, is inversely proportional to the electrochemical rate 
constant of the determining step. In the case of cathodi- 
cally treated platinum, the value of A0 is three-times larger 
than on pc platinum, indicating a large consumption of 
adsorbed species during the process. The values of A@ on 
the metal modified platinum surfaces are lower except for 
platinum/silver where the methanol adsorption is strongly 
preferred. 
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(5) The value of the electrochemical rate-determining step 
was determined in all cases. The largest value was 
observed for platinum/ruthenium/osmium surface, i.e. 
kox,4 = 1073 s~!, followed by platinum/ruthenium and plat- 
inum/osmium. The situation for platinum/silver is the worst 
with kox,4=5 x 1079s-!. These results seem to be in 
contradiction with the large initial current obtained for 
cathodically treated and bare platinum. However, these sur- 
faces are poisoned after long exposure times to methanol 
oxidation and the initial large currents get depressed. 
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Appendix A. Ricatti's equation 


The non-linear differential equation with quadratic terms 


S dh RG)y? 
dx (x) + QG)y + R(x)y 
is called Ricatti's equation. If yı is a known particular solution 


of (A1), then the substitutions of 


(Al) 


4 d dy dy du (A2) 
= and — = — + — 
da dx ^ dx | dx 
in (A1) leads to the following differential equation for u: 
du 2 
— — (Q + 2yi R)u = Ru (A3) 
dx 
Since (A3) is a Bernoulli equation with n 2 2, 
dy " 
a, ^ POY = foy, (A4) 
x 
it in turn be reduced to the linear equation 
dw 
—+(Q+2y,R)w=-—-R (A5) 


dx 


by the substitution w=u7!. The situation is that in many cases 
a solution of a Ricatti equation cannot be expressed in terms of 
elementary functions. 

For example, if we have the following differential equation 
d 
a s 2 —2xy + y 
dx 

Itcan be easily verified that a particular solution of this equa- 
tion is y; 22x. From (A1) we make the identifications P(x) 22, 
Q(x) = —2x, and R(x) = 1 and then solve the linear Eq. (6): 


(A6) 


dw dw 
— +(—2x + 4x)w = -1 or —+2xw=-1 (A7) 
dx dx 

The integrating factor for the last equation is e*, so 
d z m 
—(e =ë A8 
PA w) (A8) 


Now the integral 


x 
2 
f e! dt 
X0 


cannot be expressed in terms of elementary functions. Thus, we 
write instead: 


2 E 1 E M 
e w= -f e dt +c or e? (=) = -f e dt +c 
Xo u Xo 


(A9) 


(A10) 
so that 
e? 
u = c- fe (A11) 


A solution of the equation is then y =2x + u. 


Appendix B. The integration of the surface coverage by 
methanol adsorbates 


The time dependence of the surface coverage by methanol 
adsorbates obeys: 


dð Nae CiCx(ri — r2) exp((r1 + r2) (Bl) 
dt Tsat Fkox 4 (Cy exp(rit) +C2 exp(r21))” 


We can integrate it as follows: 


ri dis C1 C€5(1 — r2 n exp((r1 + r2)0) " 
6; kox,4 o (C1 exp(rit) + Co exp(rat))” 
(B 


or from 


0 
| dé = 
6; 


exp(rif) 
(C2 + Cy exp(ri — r)t? 


CiCo(r1 — 72)" T 
0 


Kox,4 


By defining 
C1 Cx(r} — 12)" 

ja e 2(r1 — r2) (B4) 

Kox,4 

t 

—ro)t 
ee | expiry a 
0 (Cı Co expri — r2)t) 


We can use the Substitution Theorem for integration. If 
t-g(x) is a derivable function and has an inverse function 
x=a(t), and if f(g(x))g’(x) has a primitive function, it is pos- 
sible to say that the primitive of f(t) will be f f(t)dt = 
S FE) dx. 


In our case, we have: 


exp(rı — r2)t 
(Cy + Cz expri — r2) 


If t= g(x)= 1/(rı —r2) Inx then: xzexp(ri — 72) and 
g'(x) = l/(ri — r2)x. By substituting the new variables in the 


fs (B6) 
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primitive of f(t) 


1 1 
dt = d 
"n f (r1 — r2) / (C1 + Cox)? i 


Therefore, we can integrate now as: 


(B7) 


1 | 1 F 1 | =i |. 
P — 
(ri — r2) J (Cy + Cox Crary — r2) | (C17 Cox) | 4 
(B8) 


since the integration limits are x; 2 exp(r1 — r2)t; 21 (tı is the 
Zero initial time) and x2 = exp(ri — r2)t2 


1 
doce. 
"n : Co(r1 — r2) 


x | j : (B9) 
(Cy +C2) (Ci + Coexp(ri — r2)t) 


We can obtain now the surface coverage by the final substi- 
tution: 


C = 
6=6;— 171 — r2) 
Kox,4 
1 1 
x (B10) 
(Cy +C2) (Cy Coexp(ri — r2)t) 
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Glossary 


c: the surface concentration in (mol cm?) 

F: the Faraday's constant (C mol!) 

j(t): the time dependent current density (A cm?) 

K, Land M: parametric terms of the Ricatti equation (cm? mol! s^!) with K = 
F/o[(4 = XMaeao(H*) 7 s (4 = X)kox,2 + Xkrea i (H7) vu Krea,3(H* )], 
L = Fkox,4/0, M = F/o[Kox,3 — (4 — x)kox,1(CH30H)] 
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kox,i and kyeqj: the electrochemical rate constants for each oxidation and reduc- 
tion pathways (s7!) 

n: the number of electrons 

N4: the Avogadro’s number (molecules mol!) 

pi: the reaction order for bulk and surface concentrations of reactants 

tss: the steady state time at which the variation of the current density is zero 

V: the electrochemical rate of reaction (mol cm? s7!) 

o: the maximum surface charge density in (C cm?) 

0: surface coverage by adsorbed intermediates in methanol oxidation 


0j: the initial surface coverage by methanol adsorbates (t > 0) 

Oss: the value of surface coverage at steady state (t — oo) 

A0: the integral change in the surface coverage equal to 05, — 8—0 

A: the pre-exponential factor of the current density defined as A= 
oC, C»(ri — r2) | Kox 4. Cı and C» (dimensionless) being the pre- 
exponential factors of the linear combination in the solution of Ricatti 
expression, and rı and rz (s^!) the exponential terms in the same solution 

I and T sar: the surface excess at any time and at under saturation conditions in 
(molecules cm~?) with Psat 2 NAo/F 


